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Abstract 

The spectral properties of the Laplacian on a class of quantum graphs with random metric 
structure are studied. Namely, we consider quantum graphs spanned by the simple Z d -latticc 
with £-type boundary conditions at the vertices, and we assume that the edge lengths are 
randomly independently identically distributed. Under the assumption that the coupling 
constant at the vertices does not vanish, we show that the operator exhibits the Anderson 
localization near the spectral edges situated outside a certain forbidden set. 
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1 Introduction 



In the present paper we are discussing the Anderson localization for a special class of random 
perturbations of periodic structures. As argued in the original paper by Anderson [4], there is 
no transport in disordered media under certain conditions. This phenomenon was interpreted 
initially within the framework of spectral theory: as proved in various settings, a generic random 
perturbation of periodic operators produces a dense pure point spectrum in certain energy 
intervals; at the same time one can be interested in the dynamical localization consisting in 
uniform bounds, in both space and energy, for propagating wave packets, which implies the 
dense pure point spectrum. We refer to the very recent collection of papers [11] providing 
the state of art in the theory of random Schrodinger operators, in particular, the paper [23] 
discussing various mathematical interpretations of the Anderson localization. 
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Up to now, the most studied models of disordered media are either discrete (tight-binding) 
Hamiltonians or continuous Schrodinger operators. In the last decade, there is an increasing 
interest in the analysis of quantum Hamiltonians on so-called quantum graphs, i.e. differen- 
tial operators on singular one-dimensional spaces, see the collection of papers [6, 13, 15, 27]. A 
quantum graph is composed of one-dimensional differential operators on the edges and boundary 
conditions at the vertices describing coupling of edges. Such operators provide an effective model 
for the study of various phenomena in the condensed matter physics admitting an experimental 
verification, and there is natural question about the influence of random perturbations in such 
systems [39]. There are numerous possibilities to introduce randomness: combinatorial struc- 
ture, coefficients of differential expression, coupling, metric, etc. Being locally of one-dimensional 
nature and admitting a complex global shape, quantum graphs take an intermediate position 
between the one-dimensional and higher dimensional Schrodinger operators. It seems that the 
paper [26] considering the random necklace model was the first one discussing random inter- 
actions and the Anderson localization in the quantum graph setting. Later, these results were 
generalized for radial tree configurations [22], where Anderson localization at all energies was 
proved. Both papers used a machinery specific for one-dimensional operators. The paper [1] 
addressed the spectral analysis on quantum tree graphs with random edge lengths; it appears 
that the Anderson localization does not hold near the bottom of the spectrum at least in the 
weak disorder limit and one always has some absolutely continuous spectrum. Another impor- 
tant class of quantum graphs is given by Z d -lattices. The paper [14] studied the situation where 
each edge carries a random potential and showed the Anderson localization near the bottom of 
the spectrum. Some generalizations were then obtained in [20, 21]. The case of random coupling 
was considered by the present authors in [25]; recently we learned on an earlier paper [10] where 
some preliminary estimates for the same model were obtained. 

The present Letter is devoted to the study of quantum graphs spanned by the Z d -lattice where 
the edge lengths are random independent identically distributed variables. We consider the 
free Laplacian on each edge with 5-type boundary conditions and show, under certain technical 
assumptions, that the operator exhibits the Anderson localization at the bottom of the spectrum, 
i.e. that the bottom of the spectrum is pure point with exponentially decaying eigenfunctions. 

There are two basic methods of proving localization for random operators: the multiscale 
analysis going back to Frohlich and Spencer [16] and the Aizenman-Molchanov method [3]. The 
Aizenman-Molchanov method gives explicit and efficient criteria for localization in terms of the 
Green function but only works under special assumptions on the way the randomness enters 
the problem (we used this method for the study of the random coupling model in [25]), which 
does not hold in the situation we are studying. On the other hand, the multiscale analysis is 
a rather universal tool which can handle very abstract situations [38]. It is a certain iterative 
procedure which works as far as some input data are available (see the paper [14] discussing 
the multiscale analysis for quantum graphs). Below we are concentrating on obtaining the the 
most important necessary ingredients, more precisely, the Wegner estimate and the initial scale 
estimate. We are mostly interested in the spectral interpretation of the Anderson localization, 
and our results, being combined with the multiscale analysis, prove the presence of the dense 
pure point spectrum in respective energy ranges. Nevertheless, they also can be used for the 
study of the dynamical localization; we refer to [23] for details. 

During the revision of the Letter the preprint [30] appeared, which studies a similar model 
(but with a different problem setting) and contains an alternative proof of the Wegner estimate 
for the zero coupling constants. 
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2 Random length model on a quantum graph lattice 



We recall here some basic constructions for quantum graphs. For the general theory see e.g. the 
reviews [18, 28, 29] and the collections of papers [6, 13, 15, 27]. 

Let r = (V, £) be a countable directed graph with V and £ being the sets of vertices and 
edges, respectively. For an edge e G £, we denote by ie its initial vertex and by re its terminal 
vertex. For e G £ and v G V we write f~eore~t>ifuG {ie, re}. The degree of a vertex is 
the number degf := #{e G £ : e ~ v}. 

For < £ m i n < l max < oo consider a function I : £ — > [7 m in^max]- Sometimes we will write 
l e instead of /(e); this number will be interpreted as the length of the edge e. Replacing each 
edge e by a copy of the segment [0, Z e ] in such a way that corresponds to ie and l e to re, 
one obtains a so-called metric graph. Our aim is to study a special-type differential operator on 
such a structure. 

In the space TL := (B e e£^ e > TL e = L 2 [0, Z e ] consider an operator H acting as / =: (f e ) i— > 
(— f") =: Hf on the domain consisting of the functions / G © eg £ H 2 [0, l e ] satisfying the Kirch- 
hoff boundary conditions, i.e. for any v G V one has 

fM = fb(0) =: f(v), re = ib = v, (1) 

and 

f(v)=af(v), /»:= £ £(0)- £ (2) 

e:te=i) e:re=!) 

where a is a real number, the so-called coupling constant (for simplicity, we assume that the 
coupling constants are the same for all vertices, which is sufficient for our purposes). It is known 
that the operator thus obtained is self-adjoint [28]. We denote this operator by H(T,l,a). 

We are going to study a special case of underlying combinatorial configuration, namely a 
periodic lattice with random edge lengths. Let V = Z d , d > 1, and hj, j = 1, . . . } d, be the 
canonical basis of Set 

£ d := {(m,m + hj), m£Z d , j = l,...,d}, 

where for an edge e = (v,v f ) G v,v' G one has ie := v, re := v' . For this graph 
F d := (Z rf , £j) consider the operator H (/, a) := H(T d , I, a). In the present paper, we study some 
spectral properties of such operators under the assumption that l e are random independent 
identically distributed variables. 

Namely, on (S7,P), a probability space, let ^ be a family of independent identically 

distributed (i.i.d.) random variables whose common distribution has a Lipschitz continuous 
density p with support [l m in,'max]- By a random Hamiltonian on the quantum graph, we mean 
the family of operators H u (a) := H(l u , a). 

For n G Z d consider the shifts r n acting on the set of edges, T n (m,m f ) = (m + n,m' + n) 
and the operators H(l^,a) with Z^(e) := /^(r n e). Clearly, H(l^,a) is unitarily equivalent to 
H^^a) for any n as H(l^,a)U n = U n H(l^,a), U n (f e ) = (f Tn e)- Li terms of the theory of 
random operators, the shifts r n form a measure preserving ergodic family on which allows 
one to use the standard results from the theory of random operators [9, 37]. In particular, one 
obtains the non-randomness of the spectrum and the spectral components: there exist closed 
subsets Xj = Sj(a) C M and a subset Q 1 C 17 with P(f2') = 1 such that spec^ H w {a) = Sj, 
3 ^ {pp> a c,sc}, for any lo G 0,' . We recall that the pure point spectrum spec pp i? is the 
closure of the set of the eigenvalues of H). Let £(a) = S pp U S ac U S sc be the almost sure 
spectrum of H w {a). The aim of the paper is to show that there are intervals / C S such that 
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S ac n I = S sc n / = 0; this means that the interval / is densely filled with eigenvalues of H u and 
there is no continuous spectrum in I for almost all u, or, in other words, the spectrum of H u in 
/ is almost surely dense pure point. 

Let us describe first the almost sure spectrum as a set. For u > and E M denote by 
the periodic Kronig- Penney operator acting in L 2 (IR) i.e. the point interaction Hamiltonian, 

d 2 

fcez 

The operator can be correctly defined through the associated sesquilinear form 

(f,P u ,pg) = (/', ff ') L2 +^/(h)# M ), f, g e H 1 (K), 

and is unitarily equivalent to the operator H(l u ,f3) for d = 1, where l u is the constant function, 
l u (e) = u for all e 6 E. It is well known [2] that 

spec P Ut/ 3 = { k 2 : 9 k > 0, cos /cu + sin kG [-1,1]}. (3) 

In particular, each Dirichlet eigenvalue (irn) 2 /u 2 , n = 1, 2, . . . , is a spectral edge, and the bands 
depend continuously on both a and u. As follows from the general theory of random operators 
[9, 37], one has 

E = £(«):= |J specif, a). (4) 

min /max] 

On the other hand, as shown in [35], one has the identity spec H{l u , a) := specP M a/ /^. Hence 
the almost sure spectrum S(a) of H u} (a) is a union of bands, and the bottom of the spectrum 
is given by 

k 2 : k e (0,7r// max ) and cosA;/ max + ^ sin/d max = 1, a > 0, 



infS(a) = ^ 0, a = 0, (5) 

a 

l min 



-k 2 : k > and cosh/c/ m ; n H — sinh/c/ m i n = 1, a < 0. 



Define the set 



2kd 

i I r7r 2 n 2 7T 2 n 2 i . , 

A: =U -'"I (6) 

nGZ max min 

The set consists of the spectra of the operator H(l u ) with Dirichlet boundary condition (which 
formally corresponds to a = oo) at each vertex when u ranges over the support of the random 
variables (7g ) eS £ and the point 0. Our main result is 

Theorem 1. (a) Let a/0 and Eq be an edge of E(a) n (0, +oo) that is not contained in A. 
Then, there exists e > such that the spectrum of in (Eq — e,Eq + e) DT, is almost surely 
dense pure point and the corresponding eigenf unctions are exponentially decaying. 
(b) There exists a > 1 and e > such that for a S (— oo, —a) U (— 1/a, 0) the spectrum of ^(a) 
in [inf S,inf T, + e) is almost surely dense pure point with exponentially decaying eigenf unctions. 

An immediate corollary of the above theorem and the equalities (4) is 

Theorem 2. There exist a > 1 and e > such that for a £ (— oo, —a) U (—1/a, 0) U (0, +oo) the 
spectrum of in [inf S, inf S + e) is almost surely dense pure point with exponentially decaying 
eigenf unctions. 
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Remark 3. An elementary analysis of the condition (4) shows that, if £ max — l m j n is sufficiently 
small and there are band edges of S in (0, +00) \ A. 

Our results do not establish the localization for the important case a = 0. In this case, 
both the Wegner estimate near inf T, = and the initial scale estimate (see theorems 4 and 5 
below) fail to hold. The operator H(l w ,Q) is analogous to the acoustic operator (see e.g. [32] 
and references therein); for this operator, at least in dimension one, it is known that localization 
does not hold in its strongest form at the bottom of the spectrum (see [19]). In the case a = 0, 
the reduced operator we use to study the random quantum graph is a discrete version of the 
acoustic operator. 

As in [25], the assumptions of theorems 1 and 2 also imply dynamical localization (see remark 
7 in [25]). 

3 Multiscale analysis and finite-volume operators 

Let A be a subset of edges from Denote Va := {^e : e G A} U {re : e G A} and consider the 
graph := (Va,A). Note that this graph has no isolated vertices. We will call the operator 
H\(l,a) := H(Y d ^,l,ot) the finite-volume Hamiltonian associated to A. For random operators 
with random length functions l u , we write H^(a) := H\(l u ,a) 

In what follows, we consider Hamiltonians associated with finite cubes A = A(n) constructed 
as follows: take r£N and denote by A(n) the set of edges e such that at least one of the vertices 
v G {ie, re} satisfies \v\ < n; for the corresponding set of vertices, we write V(n) := VA( n )- 

As mentioned previously, the Anderson localization for random operators can be established 
using a certain iterative procedure called the multiscale analysis. In order to start the multiscale 
analysis one needs to verify the validity of several conditions for a fixed interval Jcl, which 
then imply the localization in some subset of / in various settings, see e.g. [38, Section 3.2]. 

The first group of conditions uses very few information on the nature of random interactions, 
i.e. whether one has random edge length, random potential on the edges or the random coupling 
constants etc. and usually need only some uniform bounds for the random variables. These 
conditions £1X6 clS follows: 

(a) the finite-volume Hamiltonians and H 1 ^, corresponding to any two non-overlapping 
finite sets of edges A, A' C £d are independent, 

(b) the finite-volume operator obeys a Weyl estimate for the eigenvalues in I, i.e. there exists 
a constant C > such that the number of the eigenvalues of H% in / can be estimated 
from above by Cn d for all n G Z and almost all uj G £1, 

(c) there exists a geometric resolvent inequality which provides some uniform bounds for the 
resolvents of finite-volume operators in terms of the operators corresponding to smaller 
finite volumes, 

(d) a generalized spectral theorem (Schnol-type theorem). This means that the existence of 
a non-trivial solution / to H w f = Ef (i.e. —f" = Ef e on each edge and the boundary 
conditions at the vertices are satisfied) with a suitable bound an infinity implies E G 
spec (for E G /) and, moreover, the spectrum of H u in / is the closure of the values 
E £ I with the above property. 

Depending on the concrete problem, these conditions can weakened, see e.g. [23, Section 4]. 
Note that the first condition (a) is trivially satisfied in our case. The conditions (b) and (c) were 
shown in [14] for equilateral lattices, and the proof goes in our case with minor modifications 
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for any interval /. The generalized spectral theorem, the condition (d), holds in any interval I 
as well due to the results of [7] (see also [31] for generalizations). 

The second group of conditions are very sensible to the way the randomness enters the 
system. These two conditions are 

(e) The Wegner estimate showing that the probability for to have an eigenvalue in the 
e-neighborhood of E £ / can be globally bounded by Ce a \ A\ bd with some a,b > 1. 

(f) initial scale estimate showing that that the probability for H% to have an eigenvalue in the 
|A|~ a -neighborhood of some Eq £ / can be bounded by |A|~ 6 with some suitable a, b > 
(which depend on the dimension and other parameters). 

When all the above conditions are satisfied, the multiscale analysis shows Anderson localization 
in a certain interval around the energy Eq taken from the condition (f). 

Hence, in the present Letter, we are interested in the Wegner estimate and the initial scale 
estimates for our model, that is, (e) and (f). They imply theorem 1 by the multiscale analysis. 

Theorem 4 (The Wegner estimate), (a) Let I C (0, oo) be an interval such that IDA = 0, then 
there exists a constant C = C(I) > such that for any interval J C I, and any cube A = A(n) 
there holds 

P{specH£(a)n J ^ 0} < C|A| | J|. (7) 

(b) There exists a > 1 such that for a £ (—oo, —a) U (— 1/a, 0) the Wegner estimate also holds at 
negative energies near the bottom of the spectrum, i.e. there exists an interval I with I 3 inf S(a) 
and C = C(I) > such that for any interval J C I, and any cube A = A(n) the estimate (7) 
holds. 

Theorem 5 (The initial scale estimate). Let E be a spectral edge of H u (a) which is not contained 
in the set A defined in (6), then for each £ > and (3 £ (0, 1) there exists n* = n*(£,/3) > 
such that, for n > n* , 

P{dist (spec iT£ (n) (a),E) < i/^ 1 } < n~ 5 . 

Remark 6. By (5), for a ^ the initial scale estimate is satisfied near the bottom of the 
spectrum inf E, which is outside A. 

When a = 0, independently of the random variables 1% and the set A, the constant function 
/ = 1 satisfies H%f = 0. Hence, both the Wegner estimate and the initial scale estimate fail for 
the energy E = 0, and this is the only spectral edge; in this case, the almost sure spectrum is 
the positive half-line. 

Actually, in dimension d = 1, the operator if w (0) is unitary equivalent to the free Laplacian 
and hence shows no Anderson localization (the spectrum is absolutely continuous). Hence, one 
has a certain similarity to the Schrodinger operators with random vector potentials, where only 
localization near internal spectral edges is proved so far [17]. 

We will prove both estimates, theorems 4 and 5 by exploiting a correspondence between the 
quantum graphs and discrete operators. A similar approach was used in [25] for quantum graphs 
with random coupling constants and more details on the reduction can be found there. 

Denote by -Djf the positive Laplacian with the Dirichlet boundary conditions in L 2 [0, Zg] and 
set D% :=0 eeA Z)-. Clearly, 

specD% = J specif, specif = { : k = 1,2,...). 
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For E ^ spec-D^ consider the operators M A (l^ , E) acting on £ 2 (V A ), 
M A (l»,EMv)=VE( ]T — ^^(re)+ £ — \-= ^e) 

- cat%VE<p(vj). (8) 



Here and in what follows, the continuous branch of the square root is fixed by the condition 
Sv^E > for E £ M. 

The meaning of the operators Ma(Z £J , .£7) is as follows. Consider the equation H%f = Ef for 
E ^ specD u . On each edge e £ A, f e satisfies the Dirichlet problem —f" = Ef e , / e (0) = /(te), 
fe(le) = /( re ) ( see (-*-) f° r the definition of the values f(v)). Therefore, 

/ e (t) = /(te) ; — 7=7— + /(re)- — -=— . 

sin V sin v E/ w 

Substituting this representation into the boundary conditions (2) yields M A (l w ,E)f A = af\ 
where f A = (f( v )) v& y ■ For the complete graph, A = £ d , we simply write M(l w ,E) instead 
of M A {1^ ,E). The map E i— > M A {l u ,E) is obviously analytic outside specD 1 ^. The following 
characterization of the spectrum of H^(a) shown in [34] will be the key to our analysis: 

• an energy E ^ specD^ is in the spectrum of H^{a) if and only if a £ spec Ma (l u , E) 

• for each such E, one has dimker (H%(a) — E) = dimker (Ma(/ w , E) — a). 

For infinite A, one needs to use the self-adjoint extension theory [8]; for finite A this relation 
has been known for a long time, see e.g. [5, 12]. 

We note that similar relations between quantum graphs and discrete operators exist for 
more general boundary conditions at the vertices, but the corresponding reduced operators 
M(E) become much more complicated, see [36]. 

4 Proof of theorem 4 (Wegner estimate) 

As noted previously, one has 

P{ specific*) n = F{3E£ J: a G specM A (Z w , E)}. (9) 

Note also that, for any E* £ I, one can write M A (£ W , E) = M A (r , £*) + (E- £*) Aj Ba(^, E, £*). 
Due to analyticity, one can find a constant b > such that \\B A (l u , E, E*)\\ < b for all E*, E £ I 
and A C £■ and almost all u £ Q. 

On the other hand, the condition a £ spec Ma (l^ , E) implies the existence of a vector 
ip £ £ 2 (V A ), IMI = 1, such that (Ma(/ w , E) — a) ip = 0. Let Ej be the center of J. The estimate, 
for E £ J, 

\\(M A (l»,Ej)-a)<p\\ < \\(M A (l",E)-a)<p\\ + \E-Ej\ • \\B A (F ,E,Ej)y\\ < b\J\, 
yields the inequality 

P{spec^(a)n < p{dist(specM A (r,Ej),a) < 6|J|}. (10) 

In what follows, we denote Ej simply by E to alleviate the notation. 
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For e £ £, introduce the operators Pf , P| , I e acting on £ 2 (V\): 

1 



Pif(v) 



I e f(v) 

In terms of these operators, one has 



2 (/(«0 + /(re)), ve{Le,re}, 
0, otherwise, 
1 
2 
1 
2 

otherwise, 

_\f( v ), v£{Le,re}, 
I 0, otherwise. 



(/(ie) - /(re)), v = te, 
{f{re) - f(te)), v = re, 



,0, 



(11) 



(12) 



(13) 



M A (l",E) = £ . ; (Pf - P|) - VEcotl»VEr 



eGA * sin ^ 



and 



dM A (l",E) E cos l^VE 
3^ " ~ S in 2 l^y/E 



[Pt-Pi) + 



E 



sin 2 l^VE 



Consider first the part (a) of the theorem, i.e. the case I C (0,+oo). As ||Pf — P||| = 1 and 
pep = P? for j G {1,2}, one has 

- cos l^VE (Pf - P|) + P > (1 - | cos 1%VE\)P. 
As / does not meet A, there exist constants c±, C2 > such that 



1 — | cos l e \f~E\ > c\ and 



E 



s\n z l^y/E 



> C2 for all e G £ and E £ I and a.e. u; G f2. 



Hence, 



so that 



dM A {l^,E) 



> c\CiP for all e € £ and E £ I 



or 



D A M A (l", E)>(3id with D A := £ 

eeA e 

Let Ej^(a, b) denote the spectral projection of M A {l w , P) onto the interval (a, b). There holds 
# [specM A (r,P)n (a-b\J\,a + b\J\)) 



tr E%(a - b\J\, a + b\J\) = tr 



a+6|J| 



a-6|J| 



9 AX (-oc J o](M A (r ) P J )-A)dA 



• (14) 
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On the other hand, one has 



-tr [D AX{ - oo , 0] {M A (r,E)-X)) 

= tr [5 A X(-oo,o](M A (Z w ,^ " \)D A M A (l",E)] > /3tr [d xX{ -oo, ]{M A (l" , E) - A)] . 

The last estimate is possible as both operators under the trace sign are non-negative. Hence, 



tr [dxX(-oo,o] (M A (l",E) - A)] < /T 1 tr 
where we denoted for brevity d e := 8/81%, and 

ra+b\J\ 



J2-deX(-oo,o]{M A (l u ,E)-X) 

lee£ 



tv[E%(a-b\J\,a + b\J\)] K/T 1 / ]T tr [-0 e X(-oo,o] {M A {T , E) - A)] dX. 

<*-b\J\ ee£ 



Taking the expectation, one obtains 



Etv E A (a-b\J\, a + b\J\) < P' 1 ^2 I P( l e<) dl t> / G e (£', X 1 oj)dX (15) 



a+b|J| 



a-6|J| 



whereMA, e (/ w , I, E) is the operator M A (l w ,E) with ZJf replaced by I and 



G e (E,X,u) 



^min 



pm tr [x^oo.o] (M A , e (r, Z, £) - A)] dZ. 



As the density p is Lipschitz continuous by assumption, one can integrate by parts and obtain 

. . Wmax 

G e (E,X,u) = -p(l)F e (l,E,X,u>)\~ + / P '{l)F e {l,E,X,u)dl, 

i i — 'mm / , 

" 'min 



where 



F e (l, E, A, u) := tr X (-oc,o] (M A>e (r , I, E) — X) — X (-oo,o] (M A , e (Z w U, £) - A) 



As 8 e M A (l^ ,E) is a rank-two operator, the functions F e (l,E,X,co) are uniformly bounded by 
2. Hence, the functions G e are bounded as well, say \G e \ < G for some G > 0. Plugging this 
estimate into (15), one obtains 



/a+b\J\ 
dA = C|A||J|, C - 
-i-b\J\ 



2bG 







It remains to observe that 

dist ( spec M A (f w , E) , a) < b\ J\ } < E tr £ A (a - b\ J\ , a + b\ J\ ) . 



(16) 



Now, let us now prove part (b) of theorem 4. Below we assume a < 0. Take first an arbitrary 
interval I = (—E + ,—E_) with E + > E- > 0. Note that for E < 0, it is more convenient to 
rewrite 

M A (l", E) = J2 ( . r^7= (A e - - v^cothZ^v 7 ^/^ . 



eeA 



sinhZSV— -E 
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Then, 

and one has 
with 



d e M A {l",E) 



\E\ 



I e - cosh%V^E(Pi - Pi) 



sinh 2 \%yP-E L 
F A M A (l",E) = M A {F,E) + K A (F,E) 



Fa := --jL= ^tanh/^v^ Jj, K A (l»,E) := -L= ^tanl<^E/ e . 

* eGA e * eeA 

tanh Z ■ \/~E 

Denote 7 = ttt~ > then one has K A (l u ,E) > 7 id for all A, E E I, and a.e. cj. As in 



<E4 



the case a > 0, one computes 

-tr[F A (x(-oc,o](M A (r£)-A))] 

= tr [d x X(-oc,0](M A (l»,E) - \)F A M A (l»,E)] 
= tr [5 A X(- O o,0](M A (r, J E;)-A)M A (Z w , J E;)] +tr [a A X(-oo,o](M A (Z w , - A)K A (/ W , £)] 
> tr [Aa A X(-oo,o](M A (r,^) - A)] + 7 tr [a A X(-oc,o](M A ^, £) - A)] 

= tr [(7 + X)d xXi - 00 ,o](M A (l^,E) - A)] . 

Assume that 



A + 7 > > for A e (a - a + 6|/|) = (a - - £L|,a + - (17) 



then 



/3ti[a xX( - xfi] (M A (f,E)-k)] <tr 



e£ A 



Using this inequality and (14) one can estimate 

ra+b\J\ 



Etr [E A (a-b\J\,a + b\J\)] = tr 



a-6|J 

i «*+&|J| tanh^V 3 ^ 5 



9 A X(-oo,o](M A (r,i?)-A) dA 



eeA 



^tr X( -oo,o](^A(Z a; ,^)-A) dA 



This is then estimated exactly as in the proof of the part (a). 

Hence we need to chose the interval / in such a way that (17) is satisfied. Let E := 

— inf spec £(a). If one has mm"^ + a > 0, then (17) holds for EL < E < E + with 



(18) 



E 

\E + — E-\ sufficiently small. On the other hand, by (5), 

a = — 1d\fE tanh mn1 ^ 
2 

Therefore, it is sufficient to find values of a for which 

m ;= tanhWg _ 2 ^ tanh /minV^ > 
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A short computation shows that the sign of f(E) coincides with the sign of cosh 2 mm ^ Ed 

which is positive for E sufficiently close to as well as for E sufficiently large, and it remains to 
note that E is a monotonous function of a due to (18). This completes the proof of theorem 4. 

5 Proof of theorem 5 (initial scale estimate) 

As in the proof of theorem 4, one can show that, for some b > 0, there holds 

P{dist(speci? A(n) (/ LJ ,a), J B) < n 13 ' 1 } < P{dist(spec Af A („)(i w , E), a) < bri^ 1 }. (19) 

Hence theorem 5 is a consequence of 

P{dist(specM A(n) (L,£),a) < bn^ 1 } < n'^. (20) 

It is well known, see e.g. Section 2.2 in [38] that, in order to prove the estimates (20), it is 
sufficient to show the Lifshitz tail behavior for the integrated density of states for the operator 
M(l u , E). Note that, if E is an edge of the almost sure spectrum of -ff^(a), then a is an edge of 
the almost sure spectrum of M(Z W , E) (see e.g. [25]). Hence, it is sufficient to study the behavior 
of the integrated density of states of M(l u ,E) at the spectral edges. 

The operator M(/ tJ ,S) is closely related to the random hopping model considered in [24]; 
below, we use the very constructions of [24] and [33] to obtain the Lifshitz tails. The integrated 
density of states in our case is defined by 

_ #{A £ spec M A(ra) (/„,,£) : A < t} 
k{t) :— lim . 

n-*oo \r\ n )\ 

This limit exists almost surely and t i— ► k(t) is non decreasing. Let [pi mm , ,Li max ] be the almost 
sure spectrum of M(Z W ,£). Then, k(t) = for t < /i min and k(t) = 1 for t > 

Mmax- Denote also 

b := sup 



E 



sval\E 

By well-known arguments, see e.g. [38, Section 2.2], in order to prove (20) it is sufficient to show 

log|log[l-A;(// max -e)]| d , ,. log I log fc(pt min + e)| d 

lim ■ < and lim ■ < . (21) 

e~*o+ loge 2 e^o+ loge 2 

For n £ N define M£(E) := M(l%,E) where l%(e) = l%(e + 7) for 7 G (2n + l)Z d . By the 
Floquet-Bloch theory, the operator M^{E) admits a density of states, fe^, satisfying 

k n( E ) = TTTw ft 1 #spec M%(E,9)n (-00, E)d0 

v L 2n + l ' 2n+lJ 

where M£(E,9) differs from M^(E) only by an operator of rank at most Cn^ 1 with C > 
independent of n. As suggested in [24], in order to obtain (21), it is sufficient to show the 
analogous estimates with k(E) replaced by K(k%(E)) uniformly in n for sufficiently large n. 
Then, as noted in [33] and applied in [24], the latter asymptotics can obtained directly from the 
following local energy estimate which has been proved in [24, Lemma 2.1]. Let a E (0,6). For 
ip £ l 2 (Z d ) one has 
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where Hq is the free Laplace operator in (Z ), 



Hou(n) = ^2 (u(n) — u(m)), 



n:\m—n\=l 

and the potential W% is given by 



= E K-z^w) + E v^cot^v^ 



^sin ^v 7 ^ 

with 

/9(i) := 



\t\, \t\ > a, 
a, otherwise. 



Then, as in [2 .], following the computations done in [33], one proves Lifshitz tails for M(/ W ,i?) 
near /i max or /U m i n . This completes the proof of theorem 5. □ 
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